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O . Abstract 

m 

^ . We investigate the J/tp + c + c photoproduction in e"*" e~ coUision at the 

i-C I LEP II energy. The physical motivations for this study are: 1) such process 

was not considered in previous investigations of J/ip photoproduction in 
Oh. e~ interaction, and we show in this work that it is worthwhile to do so 

r-| ! in order to make sound predictions for experimental comparison; 2) from 

recent Belle experiment results, the process with same final states at the 
- B factory has a theoretically yet unexplainable large fraction; hence it is 

rS I interesting to see what may happen at other colliders; 3) The process can 

I be measured with high accuracy at the planed Linear Colliders(LCs); 4) 

it is necessary to take this process into consideration in the aim of elu- 
cidating the quarkonium production mechanism, especially in testing the 
universality of NRQCD nonperturbative matrix elements. We find that the 
concerned process is really important at the LEP experiment energy; within 
the theoretical uncertainties, it is of similar magnitude as the other color- 
singlet processes when transverse momentum px > 1 GeV. Nevertheless, to 
explain the recent DELPHI experimental result color-octet mechanism is 
still necessary, but with a shrunk contribution from previous analysis. And, 
it is found that J/ip + c + c photoproduction process can not be mimicked 
by the simple fragmentation scheme. 



1 Introduction 



Quarkonium physics is still an interesting research topic while the first quarkonium 
state, J/V', has been discovered for about thirty years. Due to its approximately non- 
relativistic nature, the description of the heavy quark and anti-quark system is one of 
the simplest applications of QCD. The high precise experimental results of quarkonium 
leptonic decays render heavy quarkonium to play a crucial role in investigating various 
phenomena, such as measuring the parton distribution in hadron-hadron collision, de- 
tecting the Quark- Gluon- Plasma signal and even new physics, etc. On the other hand, 
the interplay of perturbative and non-perturbative Quantum Chromodynamics(QCD) 
happens in the quarkonium production and decays, which can therefore stand as probes 
in investigating the non-perturbative nature of QCD. 

Quarkonium physics has experienced dramatic advances in recent years, among 
them the current focus in the field has been on the colour-octet mechanism(COM) [1] 
which was triggered by the high-p^- J/ip surplus production discovered by the CDF 
collaboration at the Tevatron in 1992 [2, 3, 4]. The colour-octet scenario was pro- 
posed based on a novel effective theory, the non-relativistic QCD(NRQCD) [5]. Hav- 
ing achieved the first-step success in explaining the CDF data, COM imposes as well 
a strong impact onto almost every aspects of quarkonium physics, and various efforts 
have been made to confirm this mechanism, or to fix the magnitudes of the universal 
NRQCD matrix elements. Although the theoretical framework seems to show qualita- 
tive agreements with experimental data, there are certain difficulties in the quantitative 
estimate of the colour-octet contribution [6], in particular, in J/ip and ip' photopro- 
duction at HERA [7, 8, 9, 10, 12], J/iplip') polarization in large transverse momentum 
production at the Fermilab Tevatron [13, 14, 15], and more recently in B-factories. 

It is widely expected that the B factories would provide clearer information of the 
quarkonium production. The B-factory experiments recently reported their measure- 
ments on the prompt charmonium production at e'^e~ colliders at y/s — 10.6 GeV 
[16, 17, 18]. To one's surprise, both their inclusive and exclusive measurements have 
big discrepancies with theoretical calculations [19, 20, 21, 22, 23, 24]. Among the 
puzzling features of the B-factory data, in particular, the total cross section of the ex- 
clusive + e~ J/i/j + rjc process is found to be about an order larger than theoretical 
predictions [22, 23, 24]. That is [18]: 

(7(e+ + e- ^ J/^ + ric) x Bir/c ^> 4charged) = {OmStoZl ± 0.009)pb , (1) 
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As well, the Belle collaboration [18] found a large cross section for J /if) inclusive pro- 
duction along with an open charm pair, the same final states process as what we are 
going to discuss, 

rri p'^ p~ — ^ 1 1 ill -\- n n\ 

= 0.591°:!^ ±0.12, (2) 



(T(e+ e- ^ J/V' + X) 
which is far more than theoretical expectations [19, 20, 21]. The new B-factory data, in 
some sense, pose a new "crisis" in the study of qaurkonium physics. Therefore, to reveal 
the problems lying behind the prevailing quarkonium production mechanisms (models) 
is currently an urgent task, and possibly has a long way to go. Nevertheless, the 
BELLE "puzzle" docs not really mean the failure of QCD based quarkonium production 
mechanisms, like NRQCD and CS model. The " cirsis" may stem from the unexplored 
higher order contributions, for instance, and other yet unknown reasons within the 
framework of NRQCD. 

The final establishment of NRQCD factorization as the correct theory of quarko- 
nium production and decays still needs more tests. The universality of NRQCD matrix 
elements is one of the critical points to be verified. People have tried many ways to 
discover the universahty of COM at different coUiders, as far there is still no decisive 
result either approval or disapproval of it. 

Quarkonium photoproduction in e"*" e~ collisions were investigated by several groups 
[25]. And, in very recently, based on leading order perturbative QCD analyses, Klasen 
et al. [27] find that the new DELPHI [28] data evidently favor the NRQCD formalism 
for J/ijj production, but rather the conventional colour-singlet (CS) model [29], which 
is quite encouraging. Considering the data accumulated at all LEP four detectors at 
CERN may still tell us more about the quarkonium production in the future, we re- 
alize it is meaningful to investigate the quarkonium photoproduction in more details. 
And, we find although superficially the 7 + 7 ^ J/^/' + c + c process stands as a sub- 
leading order process(in the sense of strong coupling in comparison with the resolved 
photon processes), its contribution to quarkonium production does not really neces- 
sary to be minor to other processes within the CS prescription. In the direct photon 
production, the concerned process here is obviously the leading-contribution process, 
since in experiment the J/ip + ^ final state process is suppressed. In case of resolved 
photon production, reference [27] finds out that single resolved photon processes give 
the dominant contribution. The resolved processes in general are suppressed by the 
parton distribution probability, but may get compensation from the order of coupling 
constant (s). To find whether direct or resolved proceses dominate in the J/-^ photo- 
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production at LEP, one needs to do a concrete calculation. 

The reminder of the paper is organized as follows. In Sec. 2, we will give a de- 
scription of our calculation procedure. In Sec. 3, our numerical results are presented, 
where the theoretical calculation of J /ip + c + c photoproduction at e~ collider is 
confronted with recent experimental result at LEP. Finally, we give our summary and 
conclusions. 

2 Physics Motivation and Formalism 

As explained in the introduction, we are going to supplement the process 7 + 7 ^ 
J/V' + c + c in 6+ e~ scattering to the analyses of J/'^ inclusive production at LEP IL 
Here, the colliding photons can participate in the hard interaction either directly, or 
resolvedly through their hadronic components. As realized, both direct and resolved 
processes can be in the same order within our interested energy distribution scope [27] . 
In this sense, since our concerned process is perturbatively at sub-leading order relative 
to the resolved processes considered in [27], it looks to be negligible at first glance. 
However, the complexity of this process makes the order analyses not transparent. 
And further, there is no obvious reason to disregard this process while taking color- 
octet processes into consideration. Therefore, to make an overall estimation of J/tp 
photoproduction at LEP and draw conclusions without considering the new proposed 
process may run some risk, as shown in next section. Of course, in case of having enough 
events, as inspired by B factory experiments, one may expect that the concerned process 
can be well distinguished from other processes. In addition, since the J/%l)-\-c -\- c final 
states with a relatively large invariant mass, it is easy to imagine that the resolved- 
photon contribution for this process will be less important than the direct one, which 
is confirmed by our numerical evaluation. We find in explicit calculation that the 
resolved-photon one are really negligibly small. 

On the other hand, it is easy to attribute the 7-I-7— ^J/'^-l-c -|- c process to 
a simple fragmentation representation approximately, where charmonium is produced 
via the charm quark fragmentation. It is worthwhile to mention that the situation 
here is different from and more complicated than quarkonium production in, e.g., 
decays, where the fragmentation mechanism works pretty well and the calculation can 
be greatly simplified by taking the fragmentation limit. Here, some " nonfragmenting" 
graphs are not negligible. Explicit numerical result given in next section support this 
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argument. 

Generally speaking, in photon-photon collision the interacting photons can either 
originate from the bremsstrahlung of high energy electron-position collision, beam- 
stralung, or, theoretically, be obtained by the Compton back-scattering of laser light 
off the linac electron beams and realize the photon-photon collision at a linear collider 
with approximately the same luminosity as that of the e+ beams. Here, in this work 
we will refrain our study on the first case and confront our result with the experimental 
data analyzed recently by DELPHI Collaboration of LEP II experiment at CERN. 

The source of photons from electron-position bremsstrahlung can be well formulated 
in Weizacker- Williams approximation (WWa) [30]: 



27r 



— oy-^maxi "^mmJ ' e I /02 



(3) 



where Q^^„ = m^^Vll - and g^„^ = [E ef{\ - x) + g^^„ with x = E^/E^, 9 the 
experimental angular cut in order to ensure the photon to be real, and E = Eg = y/s/2. 



Figure 1: Half of the Feynman diagrams of discussed J/ip producing subprocess, 7 -|- 
7 — > J/ip + c + c. The missing part of diagrams are the charge conjugation of the shown 
ones, and can be simply obtained by flipping the fermion flow directions. 

Our concerned process involves with twenty Feynman diagrams. Ten are shown 
in Figure 1 and the left ten are just the charge conjugation of them. It is evident 
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that in this process the quarkonium can be formed in CS configuration, which may 
be formulated coincidentally in both CS model(CSM) and NRQCD description at the 
leading order. 

The calculation of prompt J/ip producing rate will be carried out by standard 
procedure with normalization of spin project operators for the quarkonium production 
taken as: 

P _ P 11 

Vs,sAP-,q)^ Y,v{--q]S2)u{- + q]Si) < -, Si, -, S2\S, S, > , (4) 

where P and S, are respectively the quarkonium four-momentum, its spin and the 
z component of the spin; q is the relative momentum of the heavy quarks; and Si, S2 
represent their spins. In the non-relativistic limit, for S-wave states, to leading order 
the covariant forms of the projection operators are very simple: 

^o,o(P;0) = ^75(/' + M) , (5) 

Vi,sAP;0)^^f{P,S,){r + M) , (6) 

respectively, for the pseudoscalar and the vector quarkonium states. Here e^{P,Sz) 
denotes the polarization vector of the spin-1 quarkonium state, and M = 2mc is the 
mass of it. Here, rric is the charm quark mass. Projectors (5) and (6) map a QQ pair 
into the S-wave states. In our study we also repeat previous calculations, where the 
J/ip prompt production was considered; that is, the J/ip coming from higher excited 
states feeddown is taken into consideration. For P-wave states production, to leading 
order one needs to expand the relative momentum of heavy quarks to first order. Of 
the spin projectors they are 

KAP-^ 0) = — 1-[7-75(/> + M) + Tsl/' + M)7«] , (7) 

n^5.(i^;0) = ^^[7" f{r + M)+ ^(f + M)7l , (8) 

respectively. 

With the above spin projectors, amplitudes for 7 + 7 — J/'^ + c + c process 
can be obtained, which are presented in the Appendix for reference and comparison. 
Nevertheless, the matrix element squared is too lengthy to be shown here. All the cal- 
culation is evaluated by using the automatic processing code for computers, the FDC 
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[31]. Interested readers, who want to have the lengthy expressions and the correspond- 
ing Fortran program, are encouraged either to download directly from the web site, or 
write us. 

3 Numerical Results 

As stated in above, we perform the calculation of Feynman diagram algebra by comput- 
erized program, in which the spin projectors method was built in and is more suitable 
for evaluating complicated processes. The Feynman diagram. Analytic formulas and 
it's Fortran source are generated by the FDC This program was employed in past in 
calculating the J / ib electromagnetic production at electron-positron colliders [32] , and 
in many other applications. In order to further secure the applicability of this program, 
in preparing this work we repeat several other independent processes and compared 
with the results given in literature. The numerical calculation is performed in batch 
by the Monta Carlo subroutine encoded in FDC, too. 

The overall differential cross section of the J/ t/j photoproduction can be obtained by 
the double convolution of parton-parton (photon- photon) to J/ ip cross section with the 
parton distribution functions in photon and the photon distribution densities, schemat- 
ically like 

Here, the f-y{x) represents the photon density in e"*" e~ collisions or at photon colliders, 
and fi/j{x) — 7, g, u, d, s) denotes the GRS photon PDFs [33]. For direct photon- 
photon interaction, it is obvious that the f^/^{x) will be the Delta function. 

In doing numerical calculation, the general parameters are taken as: a = 1/137.065; 
< >= 1.4 GeV^; rric = 1.5±0.1 GeV; Ag^^ = 174 MeV [26] and the strong 

coupling is running with renormalization scale, n = rriT- Here, mr = + p\ is the 
normally defined transverse mass of J/'^. The factor {J/i/j + i/j') — 1.278 is used to 
include ■0' contribution. With taking the non-relativistic limit the J/ip mass is taken 
to be two times of charm quark mass, and the open charm pair is of the same mass 
as charm quark in J/ip, otherwise the gauge invariance will be broken. To keep the 
consistency with other analyses, our choice of parameters is in accordance with what 
taken in ref. [27]. For details of the choice making, e.g, the magnitudes of color-octet 
matrix elements, readers are recommended to refer to the CTEQ5L fit used by [27]. 
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Figure 2: The transverse momentum distribution of J/ip photoproduction in LEP II 
experiment. The results of 7 + 7 — > J/'0 + c + c process are confronted with the central 
values in previous study in Rcf. [27] and recent DELPHI experimental result [28]. 
The upper bound of the shaded band is obtained at renormalization scale r=0.5(Mr) 
and mc= 1.4, and the lower one at r=2 and mc= 1.6 

In figure 2 we present our process versus the other theoretical predictions [27] and 
recent DELPHI experimental result [28], where in avoiding the non- photoproduction 
process the invariant mass of 7 7 system is hmited to 1^ < 35 GeV as performed in 
experiment. The maximum angle ensuring the photon to be real in eq.(2), the 6'max, is 
taken to be 32 mrad. As shown in the figure, the previously considered leading order 
CS processes contribute less than the process we are considering here as transverse mo- 
mentum being larger than 1 GeV. When transverse momentum is small we know the 
diffractive interaction process would take the leading contribution for J / production. 
^Prom the figure, the discrepancy between experimental data and colour-singlet calcu- 
lations is shrunk after including the 7 + 7^ J/^ + c + c process, although still the CS 
contributions fall below the data even with the optimal choice of the errors. In drawing 
figure 2, FDC was used to re-calculated the former NRQCD and GSM processes done 
in [27] and got an agreement with their results numerically. 

Appearing in the figure, the results from previous study are only taken by the central 
values. We notice that there are big uncertainties remaining in both previous analyses 
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and our calculation shown as shaded band. The theoretical errors come mainly from 
the influence of scale dependence (// = O.BrriT, rriT, 2mT), the non-perturbative matrix 
element uncertainty, and the variation of charm quark mass {rric = 1.5±0.1Gey). The 
strong scale and mass dependence imply that the higher order relativistic and radiative 
corrections would be large. 

As shown in the figure 2, the 7 + 7 — >■ J/ip + c + c process can not be reproduced by 
the fragmentation mechanism [11] (simply times 2.4 x 10^'^ to 7 + 7 ^ c + c), even at 
the high px side. This finding means that the non-fragmentation-like graphs existing 
in this process arc not negligible in high energy and high pj- limit. This character hints 
that the Be photoproduction, hadroproduction as well, with the similar topology of 
Feynman graphs can not be simply replaced by fragmentation mechanism. 

In figure 3 we present the invariant mass, angular, rapidity and pseudo-rapidity 
distributions of7-|-7 J/ip + c + c process, respectively. The invariant mass of 
colliding photons starts at 6 GeV as physical requirement and ends up at 35 GeV as 
imposed in DELPHI experiment. The (pseudo) rapidity varies from -2 to 2 as performed 
in experiment as well. 



Figure 3: From the upper left to lower right, it is the plot of mass distribu- 
tion, da/dM^^; angular distribution, da /d cos 9; rapidity distribution, da/dy; pseudo- 
rapidity distribution, da/drj. 

To show the influence of the renormalization scale and charm quark mass, in Table 
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I we present the relation of total cross section relying on them. Here, the r means the 
fraction of renormalization scale on charmonium transverse mass. That is: /i — r ttit, 
where pr is the J/ip transverse momentum. The total cross sections are obtained under 
the prerequisite of: ^f{10) > pr > 1.0 GeV, W < 35 GeV, and |r7| < 2, where W and 
77 are final states invariant mass and J/ip psudo-rapidity, respectively. 



TABLE I. Renormalization scale and charm quark mass dependence of the total cross- 
section. 





me = 1.4 GeV 


me = 1.5 GeV 


me = 1.6 GeV 


r = 0.5 


0.82 pb 


0.54 pb 


0.37 pb 


r = 1 


0.47 pb 


0.31 pb 


0.21 pb 


r = 2 


0.30 pb 


0.20 pb 


0.14 pb 



It is evident that both renormalization scale and charm quark mass induce large un- 
certainties on total cross section. Among them, the scale dependence is basically more 
important. Under the same physical cut and parameter input we find the new process 
at LC(s), e.g. TESLA, gives a total cross section of 0.41 pb at the central values of 
scale and charm quark mass, which is larger than the total cross section at LEP as 
shown in Table 1. This means that the new process can surely be measured with high 
accuracy at future LCs. 

4 Discussion and Conclusions 

In this work we proceed the calculation of 7 + 7— ^J/V' + c + c process of the J/ ip 
photoproduction in e'^ e~ interactions in LEP II experiment, which was not taken 
into consideration in previous analyses. We find this process is quite large in the 
photon-photon coUision at LEP. The importance of including this process hes in two 
aspects: 1) it is the dominant precess among all the CS subprocesses with transverse 
momentum greater than 1 GeV. Since in the low pr region the diffractive interaction 
will be overwhelmingly large, to have a clear cut of it one needs to focus on the large 
Pt area; 2) considering the large uncertainties still remain in the Color-Octet matrix 
elements, in goal of either to fix these uncertainties, or to test the accuracy of QCD 
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perturbative calculations, including this new process is very necessary. 

The new process we considered is unique relative to other processes of J produc- 
tion in e"*" e~ scattering. It should be obvious to obtain gauge invariance for the total 
amplitude. We performed the check and found there is gauge invariance. It is also 
checked for other processes in the aim of convincing us more of our results. In practice, 
we performed the gauge invariance check in two different ways. One at the amplitude 
level, we replace the photon polarization vector(s) by the corresponding momentum, 
then reallocate the independent terms and then numerically calculate the amplitude 
square; in another way, we do the replacement at matrix clement square level and then 
evaluate it. In both cases the large cancellation happens, and up to the precision limit 
of Fortran program they get to be zero. 

We found that the uncertainty induced by scale variation is quite large, which 
means the higher order corrections could be big. The uncertainties remaining in the 
quarkonium non-perturbative matrix elements and charm quark mass are also the 
sources of theoretical predication errors. 

We also compared the pure fragmentation result with our full calculation, and find 
that the situation of quarkonium production here differs from what in decays, where 
the fragmentation scheme can almost reproduce the full Feynman diagram calculation 
results. It is noticed that after considering our proposed process, the preliminary 
DELPHI data are still not explainable by CSM alone, and colour-octet scenario is still 
necessary. Nevertheless, since large uncertainties remains in both CS and NRQCD 
analyses, quantitative conclusions for the universality of colour-octet matrix elements 
still hard to get. In our opinion, to have a full NLO calculation of prompt J/%1) pro- 
duction would be critical on this point and beyond. Although the present DELPHI 
data are just marginal in observing the 7-I-7— > J / %l) -\- c -\- c signal, it would be still 
interesting for experimenters to see whether the accumulated LEP data are enough 
or not to find it. Anyhow, at LCs this new process should be observed with a high 
precision. It is also noticed that in previous NLO calculation of J/%1) photoproduction 
at HERA [12, 34], the similar same order process 51 + 7 —> J/V'('0') +c-|-c was missing. 
Whereas, naive estimation tells us the missing part should be not so important as the 
case in photon-photon interactions. Detailed investigation on this will be presented 
elsewhere. 
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Appendix 

The matrix element of process 7 + 7 J/ip + c + c. Here, ei, £2 £3 are the polariza- 
tions of initial photons and J/ip, respectively; and Pi{i = 1, 2, 3) are the corresponding 
momenta of them. 

M = c M(p4)(ciP3e3Piei + C2P3e3 + C3e2PiP3^\ + CAe2Pie^ei + c^e2P2.ezPiei + ceeVsesPi 
+C7e2P3e3ei + Cse2ei + cioe2P3ei + £1162^363 + 01262^3 + £13626361 + 0146263 + ci562Pi6'~i 

+01662^163 + £1762 + CisCsPie'i + Ci9636\ + 02063 + C21 + C22Pltl + C23P3j5l6'l + C24|>36'l 
+C25P3 + C26ei + C27P1 + C28P3P1 + £2963^1 + C^oPlPziz + C3162P1P3 + C3262Pi)w(p5) 

(A-1) 

where 

c = E,.i,2o y, (A-2) 



yi = {4:X3Xio{2ml - Xg - 2xio + 2^15))"^ I/2 = (-2x3(2m^ + xi5){Aml + 2x15))"^ 
y3 = {2x3Xs{2ml + Xig))"^ y^ = {2x3Xs{2ml - X2 - 2x3 + Xu))~'^ 

1/5 = (-2a;4(2m2 + a;i4)(4m2 + 20:14))-^ ye = [2x^x^{2ml + 0:14))-^ 

y-j = (2x4X8(2771^ - X2 - 2x4 + Xi5))~^ y^ = (-X8(2m^ + Xi5)(2m^ + xu + X15 + 2x19))"^ 
y9 = (-X8(2m^ + X14) (2xi - X2 - a^s))"^ Z/io = (-X2(2m^ + xu) (2xi - X2 - x^))'^ 

(A-3) 

X\=V\- P2, X2=Pl- Ps, X3=pi- P4, X4 = pi ■ P5, X5 = pi ■ 61, Xq = Pl ■ 62, 
Xj =Pl- 63, X8=P2- P3, Xg=P2- P4, Xiq = P2 ■ P5, Xn = P2 ■ 61, X12 = P2 " 62, 
2^13 = P2 ■ 63, Xu = P3 • P4, 2^15 = PS ' PS, ^^le = P3 " 61, X17 = P3 • 62, X18 = Ps ■ 63, 
2^19 = P4 ■ Pb, X20 = P4 ■ 61, X21 = P4 ■ 62, X22 = P4 • 63, X23 = Pb ' 61, X24 = Ps " 62, 
X25 — Pb ■ 63, X26 = 61 • 62, X27 = 61 • 63, X28 = 62 • 63, 



Cl = (4x24(4l/io + 1/3 + 3y4 + 3l/6 + 1/7 + 3|/8 + 4|/9) + 4X21 (-1/3 +1/4 + 1/6 

-1/7 - l/s) + 4xi7(2yio + 2y4 + 2^6 + 1/8 + 21/9)) 



(A-4) 
(A-5) 



C2 = (a;26(4x8(-l/8 - 1/9) + 8X4(1/10 + 1/6 + 1/7 + 21/9) + 8X3(-|/6 - 1/7 - I/9) 

+4^2(1/10 - 2?/6 - 2?/7) + 8xi(-yio + 1/6 + 1/7) + 8(2mc^i/8 - Xioi/s 

+2^151/8)) + (8X242^23(31/6 + 1/7 + 1/8) + 8X230:21(1/6 - 1/7 + l/s) + 8X242:20 ( (A-6) 

-1/3 - 3i/4 + l/s) + 8x21X20(1/3 - 1/4 + l/s) + 4x24a;i6(-l/io + 3^/8 - I/9) 

+4x2lXi6(l/l0 + 1/8 + 1/9) + 16(-Xi7X20l/4 + XnX2'iy%))) 

C3 = (4x7(?Jio - 3?/3 - ?/4 + y6-y7 + 1/9) + 4x25(-2yio + 1/3 " 1/4 - 1/6 + 1/7 " 3l/8 - 2^/9) 
+4X22(22/10 + 32/3 + 1/4 + 1/6 + 3^/7 + 1/8 + 2^9) + 4X13(7/10 - 27/3 " 2l/7 + I/9)) 

(A-7) 

C4 = (16mc^(-2y3 - 2^7 + ys) + (4x8(-yio + 2?/3 + 27/7 - yg) + 4x2(-yio 

+3?/3 + 1/4 - 1/6 + 1/7 - 1/9) + 4xi5(2yio - 1/3 + 1/4 + 1/6 - 1/7 + 37/8 + 21/9) (A-8) 
+4x14 (-27/10 - 37/3 -yi-y%- 3y7 - ys - 27/9))) 

C5 = 4mc(-yi - yio - y2 - ys - y4 - y^ - ya - yi - ys - yg) (A-9) 
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C6 = (4a;23(z/io - 4^6 - 4y7 + Sys - yg) + 4a;2o(-3yio - 4^6 - (A-IQ) 
-ys - Syg) + 4xi6(-yio - 4^6 - 4y7 + y8 - Syg)) 

C7 = (4x8(y8 + yg) + 8x4(-yio + ye + y? + Vs) + 4x3(2yio + 2ys + 2y4 

+4y6 + 4y7 + yg + 4yg) + 4^2 (4y6 + 4y7 + 3yg) + 4xi(yio - 4y6 - 4y7 (A-11) 
-Syg) + 4(-7mc^y8 + 5xioy8 - 5xi5ys - Sxigys)) 

Cg = (mc^(16x7(yio + 4y3 - yg + yg) + 96x25y8) + {^X25X2{-y3 + ?/4 - Vs) 

+8x220:2 (-yio - 3y3 -Vi-Vs- yg) + ^x^xi^iy^ - y^ + ys) + 8xrXu{ /^_-^2) 
yio + 3y3 + y4 + ys + yg) + 8(3xi3,Ti5y8 + 2a;i3a;2y3 + ^xuX25y8 

-3Xi5X22y8 - 3X25X8y8 - '^XjXsysjj) 

Clo = 8mcX7(yio + ys + ye + y? + yg) (A-13) 

cii = mc(8a;23(-y5 - Ve - Vr + ys) + 8x2o(yi + y2 + ys + y4 + ys) + 4xie(yio + ys + yg)) 

(A-14) 

C12 = (X27(4x8(-y8 - yg) + 8x4(2yio - ye - y?) + 8a;3(-yio - 3ye - 3y7 - 3yg) 
+4a;2(yio - 6y6 - 6y7 - 4y9) + 8x1 (-yio + 3y6 + 3y7 + 2y9) + 8(2mc^y8 
-alloys + a^isys)) + (8a;7a;23(-2yio + 2y6 + 4y7 - 3y8) + 8a;25a;23(y6 
-y7 + ys) + 8x2zX22{-y& - 3y7 + 3y8) + 8x7X2o(yio - 3y3 - y4 
+3y6 + 3y7 + 3yg) + 8a;25a:;2o(y3 - ^4 - ys) + 8a;22X2o(3y3 + y4 + ys) 
+4x7a;ie(-yio + 6y6 + 6y7 - ys + 3yg) + 4x25X16 (-yio - ys - yg) 
+4x22a;i6(yio + ys + yg) + 16x23a;i3(y7 - ys) - 16xi3X2oy3)) 

(A-15) 

cia = (mc(12x8(-y8 - yg) + 8x4(yio - 2y6 - 2y7 - y8 - yg) + 8x3(yi - 2yio -y^-y^ 
-3ye - 3y7 - 4yg) + 4x2(-3yio - 2y5 - Sye - 8y7 - 2y8 - lOyg) + 8x15 (4y8 + yg) 
+24xi(y6 + y7 + yg) + 8(-3xioy8 + a:;i4yg + ^^igys + Xgyg)) + 8mc^{7ys + 4yg)) 

(A-16) 

ci4 = (mc^(64x23(y7 - ys) + 16(-Xi6y8 - 4x2oy3)) + (i6x8X23(-y7 + ys) 

+8x23a:;2(yio - 3ye - 5y7 + 3y8 - yg) + 8x203:^2 (-2yio + 3y3 + y4 - 4y6 - 4y7 

-4y9) + 4x8Xi6(y8 + yg) + 8x4X16 (-2yio + ye + y?) + 8x3X16 (yio + 3y6 
+3y7 + 3yg) + 4x2Xi6(-2yio - 2y6 - 2y7 + yg - yg) + 8x230:15 (-ye + y?) 
+8x20X15 (-y3 + y4 + 2yg) + 4xi6Xi5(yio + yg + yg) + 8x232:14 (ye + 3y7 
-3yg) + 8x2oa;i4(-3y3 - y4 - ys) + 4xi6Xi4(-yio - ys - yg) + 8xieXi(yio 
-3ye - 3y7 - 2yg) + 8(xioXieys + 2x2oa:;8y3))) 

(A-17) 

Ci5 = mc(8x7(yi - 2y3 + ye - y?) + 8x25 (-yio + I/2 + ys - 1/4 - ye + I/7 - 2ys 
-yg) + 8x22(-yi + 2yio + 2y3 - y5 + 2y7 + yg + 2y9) + 8xi3(-y3 - y7)) 

(A-18) 

cie = mc(8x23(yio-3ye-3y7 + 2yg-y9) + 8x2o(-2yio-3y6-3y7 (A-19) 
-yg - 4yg) + 4xie(-3yio - 2y5 - 8ye - 8y7 + yg - 7yg)) 
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{x27{mc{Sx8{y8 + Syg) + 16a;4(-2yio + ye + y?) + 32x3(^6 + y7 + yg) 
+16a;2(-yio + Sye + 2y7 + 2y9) + 16xi5(-y8 - yg) + 16xi(yio - 2y6 - 2y7 
-yg) + 16(a;ioy8 - .^Myg)) + 32mJ^{-y8 - 2yg)) + me (16x70:23 (ye - y? + ys) 
+16x25X23 (-ye + y7 - ys) + 16x23X22 (-ys + 2y7 - Syg) + 16x7X20 (-yi 
+yio + 2y3 + yg) + 16x250:20 (-y2 - ys + y4 + 2y8) + 16x220:20 (yi - 2y3) 
+8x7Xi6(yio + ys + yg) + 8x25a:;ie(yio + ys + yg) + 16x23Xi3(-y7 + 2y8) 
+ 16x20Xi3(y3 - yg) + 8(xi3Xi6ys - 2xi6X22ys))) 

(A-20) 

= mc(8x24(4yio + ys + 3y4 + 3ye + y7 + 2y8 + 4y9) + 8x21 (yio + 2y4 + 2y6 
-ys + yg) + 4xi7(5yio + 2y2 + 3y3 + 5y4 + Sye + 3y7 + Sys + 5yg)) 

(A-21) 

(8mc^(5xi7y8 - 4x24y8) + (16x24X2(2yio + 2y6 + y7 + ys + 2y9) 
+8x2iX2(yio + 2y6 + yg) + 4x8Xi7(-y8 - yg) + 16x4Xi7(-y7 - y8) 
+8x3Xi7(-yio - ys - y4 - ye - 3y7 - ys - yg) + 4x2Xi7(4yio + 4y6 (A-22) 
-4y7 + ys + 3yg) + 8x17X1 (ye + 3y7 + yg) + 8(-4xioXi7y8 - xi4X24ys 
+3xi5Xi7y8 - xi5X2iy8 - 2xi5X24ys + 3xi7Xi9y8 + X24a:;8y8))) 

= (a;2e(?"c(8x8(-y8 - 2yg) + 16x4(yio + ye + y? + 2yg) + 16xi5(y8 + yg) 
+16xi(-yio - yg) + 8(-2xioys + 2xi4yg + Xayio)) + 32mc^(y8 + 2yg)) 
+772^(16x24X23 (3ye + y7 - ys) + 16x24X2o(-ys - 3y4) + 16X2lX2o(-2y4 
+y8) + 8x23a;i7(5ye + 3y7 - Sys - 2yg) + 8x2oa;i7(-2y2 - 3y3 - 5y4 
-y8 - 2y9) + 8x24a:ie(-3yio - ys - 3yg) + 16x2iXie(-yio - yg) 
+4xi7Xie(-5yio - 3y8 - 9yg) + 32x2iX23ye)) 

(A-23) 

(X26 (32772^^X7 (yio - 2y6 + 2y7 - 2y9) + (16x8X7(ye - y7 + yg) + 32x25X2( 
ye + yg) + 16x22X2(-yio - 2y7) + 32x7Xi5(-ye - yg) + 16x7,Ti4(yio + 2y7) 
+16x2Xi3(-y6 + y7 - yg))) + X27(mc^(32xi7(-y8 + yg) + 64x24y8) 
+ ( 16x24X2 (-3yio - 4ye - 2y7 - ys - 3yg) + 16x2iX2(-yio - 2y6 - yg) 
+16x3Xi7(yio + ye + y7 + yg) + 8x2Xi7(-3yio - 4ye - 3yg) + 8xi7Xi5(-2y8 
+y9) + 16x17X1 (-ye - y7 - yg) + 8(2xioXi7y8 + XuX^yg + 2xi4X24y8 

+2Xi5X24ys + Xi7X8y8 - 2x24X8y8))) + X28(mc^ (32X23 (ye - y? + ys) 

+32x2o(ys - y4 + ys) + 16xie(yio + 3y8 - yg)) + (16x23X2(yio + 3y7 - ys 
+yg) + 16x20X2 (2yio - ys + ye + 3y7 + 2yg) + 16x4Xie(yio + ye + y? + yg) 
+8x2X16 (4yio + 2ye + 6y7 - ys + 4y9) + 8xieXi5(2y8 - yg) + 16x23X14 (-y? 
+y8) + 16x2oXi4(y3 + ys) + 8xieXi4(yio + ys) + 8(-2xiXi6yio 

-2xioXiey8 - 2xi5X2oy4 + 2xi5X23ye - xieXgys + 2x2oX8y4 - 2x23X8ye))) 

+ (16X7X2sXi7(-3y7 + yg) + 16x23X22Xl7(y7 - ys) + 16x7X20Xi7(-yi0 

+ys - ye - 3y7 - yg) + 16x22X2oXi7(-y3 - yg) + 16x7X24X16 (2yio + 3ye + y7 
+y8 + 2y9) + 16x7X21X16 (ye - y?) + 8x7X17X16 (yio + 2y6 - 6y7 + ys + yg) 

+8X22Xi7Xi6(-yi0 - ys - yg) + 16(xi3XieX24ys - Xi3Xi7X2oy4 
+Xi3Xi7X23y6 - Xi6X22X24y8 - XieX24X25ys + Xi7X2oX25y4 - Xi7X23X25y6))) 

(A-24) 
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C22 = (x28(16mc^(?/io - ?/3 + ?/4 + ye - yr + ?/8 + yg) + (8x8(-yio - y4 - ye - yg) 
+8x2 (-yio + ys - ye - yg) + 8x15 (yio + y4 + ye + 2y8 + yg) + 8x14 (-yg 

-y? - ys))) + (8a;7a;i7(yio - ys + ye + yg) + ^X2^xi'j{-yio - - 
-2y8 - yg) + 8a;22a;i7(y3 + y? + y8) + 8x17X13 (yio + y4 + ye + yg))) 

(A-25) 

C23 = 8a;28mc(-y2 - Vi - y? - ys) (A-26) 

C24 = (a;28(4a;8(y8 + yg) + IQx^iyj + yg) + 8x3 (yio + ya + y4 + ye + 3y7 + ys + yg) 
+4x2 (yio + 2y6 + 6y7 + 2yg) + 8x1 (-yg - 3y7 - yg) + 8(-5mc^y8 + 4xioy8 
-4xi5y8 - Sxigyg)) + (16x7X24(-2yio - 2y6 - y? - Vs - 2y9) + 8x7X2i( 
-yio - 2y6 - yg) + 4x7Xi7(-5yio - Gye - 2y7 - ys - Syg) + 8(-xi3X24y8 
+xi7X25y8 + X2iX25ys + X22X2Ay?, + 2x24a:;25ys))) 

(A-27) 

C25 = "^ca:2s(16x23(-y7 + ys) + 16x2o(y2 + y3 + ys) + 8.Ti6ys) (A-28) 

C26 = (x2s(mc(8x8(2|/8 + yg) + 32x4(y7 + y8) + 16x3(yio + y4 + ye + 2y7 + ys + yg) 
+8x2(yio + 2y6 + 6y7 + 2ys + 2yg) + 16xi(-y6 - 2y7 - yg) + 16(3xioy8 
-Xi4y8 - 4xi5y8 - Sxigys)) - 112mc^y8) + mc(16x7X24(-3yio - 3y6 - y7 
-21/8 - 3y9) + 16x7X2i(-yio - 2y6 - yg) + 8x7Xi7(-4i/io - Sye - 3y7 - ys 

-4|/g) + 16(-Xi3X24y8 + Xi'jX25y% + X2lX25y8 + X22X2m + 2X24a;25y8) ) ) 

(A-29) 

C27 = (mcX27(16x24(3yio + 3y6 + y7 + ys + 3yg) + 16x2i(yio + 2y6 + yg) 

+8x17 (4yio + Sye + 3y7 + ys + 4yg)) + mcX28( 16x23 (-ye - 2y7 + ys) 
+16x2o(-yio - ye - 2y7 - yg - yg) + 8xi6(-yio - 2^/6 - 6y7 - yg)) (A-30) 

+m,,,T26(16,T7(?/6 - ?/7 + t/g) + 16X25(-y6 + y? " yS " 2yg) + 16X22( 
yiO - y5 + 2y7) + 16Xi3(-l/7 + yg))) 

C2S = (a;27(8x24(3yio + 3y6 + y? + ys + 3yg) + 8x21 (ye - y7) + 4xi7(3yio + 4^/6 

+ys + 3yg)) + X28(8x23(-yio " ye " 3y7 + 2^8 - yg) + 8x2o(-2yio - ye - 3y7 
-ys - 2yg) + 4xi6(-3yio - 2y6 - 6y7 + ys - 3yg)) + X2e(8x7(y6 - y7 + yg) 
+8x25(-ye + yr - ys - 2yg) + 8x22(yio + ye + 3y7 + yg) + 8xi3(-2y7 + y^))) 

(A-31) 

C2g = (x2e(16mc^(-yio + 2yQ - 2y7 + ys + 3y9) + (8x8(-ye + y? - 2yg) + 16x2( 
-ye - yg) + 8x15(21/6 + ys + 3yg) + 8xi4(-yio - 2^/7))) + (8x23Xi7(y6 + 3y7 
-2ys) + 8x2oXi7(yio + ye + 3y7 + ys + yg) + 8x242^16 (-3yio - 3y6 - y7 
-ys - 3yg) + 8x2iXi6(-ye + y7) + 8xi7Xi6(-yio - ye + 3y7 - ys - yg))) 

(A-32) 

C30 = 8x26mc(-y5 - ye - y7 - yg) (A-33) 
C31 = 8x27mc(yio + ys + ye + y7 + yg) (A-34) 

C32 = (x27(16mc^(yio + 4^7 - ys + yg) + (8x2 (ye - y7) + 8x15 (-yio - ye 

+y7 - ys - yg) + 8x14(2^/10 + ye + 3y7 + 2yg) - IQxsyj)) + (8x7X16 (yio 
+2y7 + yg) + 8x25Xi6(-2|/7 + ys) + 8x22Xi6(-3yio - 2^6 - 4y7 - 3yg) 
+8xi6Xi3(yio + ye + 3y7 + yg))) 

(A-35) 
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